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0\ ' Abstract 

^ ! A set of quantitative predictions for the D-N asymmetry in the Su- 

I ; per - Kamiokande detector is presented. For these predictions, neu- 

trino events are collected in "samples" defined by the trajectories in 
^ ! the Earth of the corresponding solar neutrinos. The more important 

samples considered here are: the full sample of neutrinos detected at 
night, the sample of night neutrinos which cross the Earth core and 
^ the sample of night neutrinos which does not cross the Earth core. 

5^ I For energy integrated event rates, the D-N asymmetry for core cross- 

ing neutrinos is up to six times bigger than the D-N asymmetry for 
the full sample of night neutrinos. When the reduction in statistics is 
considered, an effective enhancement up to a factor ~ 2.3 is obtained. 
The selection of core crossing neutrinos may be relevant for the recoil 
- e~ spectra too. From these results it is reasonable to expect the 
Super - Kamiokande detector be able to constrain the small mixing 

angle solution of the solar neutrino problem (sin^ 29y < 0.01) through 
the D-N effect. Furthermore, it is proposed that the D-N effect may 
limits the total flux of neutrino produced in the Sun. At last, 
the sensitivity of these predictions to the Earth model uncertainties 
is quantitatively discussed. 



1 Introduction 



This talk reports the main results from a research conduced by the speaker 
in collaboration with Serguey T. Petcov (SISSA - Trieste, Italy, also at Insti- 
tute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 
1784 Sofia, Bulgaria) and Q.Y. Liu (SISSA - Trieste, Italy) about the Super - 
Kamiokande sensitivity to the "day-night effect" 0] . The research was aimed 
by two main purposes: to produce quantitative predictions for the day-night 
effect as seen by the Super - Kamiokande detector, to study the enhancement in 
the day-night asymmetry obtained selecting those neutrino events due to solar 
neutrinos which cross the Earth core in their way to the detector. 

The observation of the Earth or day-night (D-N) effect would be a proof of 
the validity of the MSW solution of the solar neutrino problem, since the Earth 
effect in neutrino propagation is a direct consequence of neutrino oscillations 
in matter In its simplest version the MSW mechanism involves matter - 
enhanced two-neutrino transitions of the solar z/g into an active neutrino or 
z/^, z/g ^fi{T)^ while the propagates in the Sun. In this case the ^^i{t) 
transition probability depends on two parameters: Am^ and sin^ 2^v^, where 
Am^ > is the neutrino mass squared difference and 6y is the angle charac- 
terizing the neutrino mixing in vacuum. If Am^ ^ and 6^ ^ 0, the neutrino 
propagation is sensitive to the matter distribution along the propagation path 
and the probability to detect at Earth a produced in the Sun is a function 
of the detection time t because at night the Sun is below the horizon and the 
solar neutrinos cross the Earth reaching the detector, so that further z/g ~^ ^fi{T) 
transitions occur. This can lead to a difference or "asymmetry" between the 
signals caused by the solar neutrinos in a solar neutrino detector during the 
day and during the night, i.e., to a D-N asymmetry in the signal. No other 
mechanism of depletion of the solar z/g flux proposed so far can produce such 
an effect. 

For the Bahcall and Pinsonneault (1995) solar model and assuming MSW 
i^^{t) transitions the solar neutrino data can be described for val- 
ues of the parameters Am^ and sin^ 26v belonging to two intervals [0] : 

3.6 X 10-^ eY^<Am^<9.8 x 10"^ eV^ 4.5 x lO^^ < sin^ 2^^ < 1.3 x lO^^; or 

5.7 X 10-6 eV^<Am2< 9.5 x 10"^ eV^ 0.51 < sin^ 2^v^ < 0.92. These inter- 
vals correspond to the nonadiabatic (small mixing angle) and to the adia- 
batic (large mixing angle) solutions of the solar neutrino problem. When 
theoretical uncertainties in solar neutrino flux components are considered, 
larger "conservative" regions [g, g]: 3.0 x lO"*" eV^ < Am^ < 1.2 x 10"^ eV^, 
6.6 X 10-4<sin2 2^y < 1.5 x lO^^; or 7.0 x 10"^ eV^ < Am^ < 1.6 x 10"^ eV^ 
0.30< sin2 2^^< 0.94. 



To reduce the complexity of the required calculations experimental details 
such as the background subtraction and the systematic errors were neglected, 
also continuous data taking, with no interruptions, were assumed. The only 
experimental constraints which must be considered are: the 5 MeV threshold 
in the recoil-e^ kinetic energy T^^th: and the total yearly event rate for the given 
threshold energy, which is expected to be ^ 10^ events/year which allows a 
statistical error « 1% for one year of data. In this way, predictions presented 
here are for an idealized Super - Kamiokande detector at its best performances. 
They should be considered the starting point to produce predictions which take 
in account background, systematic errors and the data taking history. Apart 
from the Bahcall and Pinsonneault (1995) solar model and the Stacey (1977) 
Earth model, other source of data are [g]. 



2 Calculating the Earth Effect 

The probability for the process Ue for a solar neutrino propagating in the 
Sun and Earth is given by [0, [T^] 

Pe(^e - ^e) = Pdl^e ^ ^e) + ^ ' ^^^^^^^ ""^^ {p^, _ sin^^^), (1) 

COS ZUv 

where /^©(i^e i^e) is the average probability of solar z/g survival in the Sun 
and Pe2 is the probability of the 1/2 transition after the have left the 

Sun. The probability Pe2 is not a constant in time because it is a function of 
the trajectory followed by the neutrinos crossing the Earth, determined by the 
instantaneous apparent position of the Sun in the sky. Given the fact that data 
are taken over a certain interval of time, the instantaneous Pe2 in eq. (j^) has to 
be replaced with its time averaged, < Pe2 >, which requires the computation of 
the solar position in the sky. 

For the scopes of D-N effect calculations it is possible to assume for the 
Earth structure a spherical symmetry, so that the electronic density is radially 
distributed. In this case the neutrino trajectory and the serie of geological 
stratifications it crosses reaching the detector are defined by the Nadir Angle 
for the Sun /i, which is the angle between the direction of the Earth center 
(Nadir) and the direction of the Sun as seen by the detector. From the detection 
time i of a neutrino event it is possible to recover h{t) allowing the selection 
of events in accord with the trajectory of the incoming neutrinos to form a 
set of five "samples" labeled: Day, Night, Mantle, Core and Deep - Core. All 
the quantities which refers to a specified sample (probabilities, spectra, etc.) 
are labeled by an index s = D, N, C, M, DC for Day, Night, Core, Mantle 



and Deep - Core samples respectively. Samples are defined as follow: the Day 
sample is formed by all neutrino events detected at day; the Core sample by 
those neutrinos which cross the Earth core {h < 33.17°); in the Mantle sample 
neutrinos does not cross the Earth core; {h > 33.17°); the Night sample is 
defined merging the Core and Mantle sample; at last in the Deep - Core sample 
neutrinos cross the inner 2/3 of core. 

The computation of Pq{V(, v^) is accomplished following the prescriptions 
derived in ref. []TT[].For a given Earth model and a fixed set of values for the 
parameters h, s\v?26v and E^j /Sm^, Pe2 is computed solving numerically the 
ordinary differential equations for the neutrino propagation in the Earth. The 
differential equation is integrated through a order Runge - Kutta algorithm 
with an adaptive step size. A set of tests assures this procedure is stable and 
accurate at the « 10"^ level. Furthermore, all over the program checks were 
made to assure consistency and proper numerical error control, based on the 
comparison of results obtained by the FORTRAN code with results obtained by 
a MATHEMATICA code. The comparison shows relative differences « 10~^ in all 
the relevant cases. In conclusion the numerical error is under control and quite 
negligible. 

Apart from the Day sample, each sample s = C, DC requires a differ- 

ent calculation of the time average probability < Pe2 >^ from the instantaneous 
probability Pe2- For a sample s = D, N, C, DC defined by the nadir angle inter- 

S S 

val [hi, /i2], and obtained over an exposure time interval [Tf, T2] (assumed here 
to last from January 1^^ 1996 to January 1'^* 1997) the averaged Pe2 probability 
is defined as: 



s ^ , , ^ s. 



where: (5 is for h outside the nadir angle interval and 1 for h inside it; R{t) 
is the Sun - Earth distance expressed in units of the mean distance (the Astro- 
nomical Unit Rq = 1.4966 x 10^ km); and 

res Jj.s ^2(^^^J V J 

s ^ s 

is the residence time, i.e., the time spent by the Sun in [h^, 7^2] • In the absence 
of neutrino oscillations the total fraction of solar neutrinos induced events in 

S ^ S 

[hi,h2] is T^gg/Ty, where Ty = {T2 — Tf) in this case is the year length. For 
the Super - Kamiokande detector the times over one year of data taking are: 
Trl/Ty = TgjTy = 0.5, r,^,/r, = 0.071, T^/Ty = 0.429, T.^f/T, = 0.0394. 



With a statistics of 10^ events this corresponds to a relative statistical error of 
about 1.4% {Night 01 Day), 3.8% (Core), 1.5% {Mantle), 5.0% {Deep - Core). 
The residence time is also relevant because it weights the contributions to 
< Pe2 >^ from each geophysical structure and from the error source. 

Many types of systematic errors can affect the predictions for the D-N asym- 
metry. Since the initial hypothesis, the most relevant error sources for these 
calculations are related to the Earth model and the apparent solar motion ||Tj. 

Earth model uncertainties introduce a relevant error source (especially for 
the core) since they affects the Earth electronic density Ue oc Y^p, where Y^, is 
the fraction of electrons per nucleon and p is the matter density (in Dr/cm^), 
the combined Ue uncertainty is « 7% -!- 10% 0. Its effect was studied repeating 
the calculations taking the p distribution [j7| but for two different hypothesis 
about in the core: = 0.467 which is representative of current knowledge 
and = 0.5 which combines all the uncertainties in Earth models. In both the 
cases it is assumed = 0.5 in the mantle |2j]. 

The apparent solar motion is quite complex and approximations are required 
to save computing time. The limit in the accuracy for the Sun's tracking is fixed 
by the angular diameter of the neutrino production region ~ 3 arcmin []T|. 
There are two kinds of approximation made which are relevant for the Earth 
effect. In the first the Earth motion around the Sun is circular (COA), in the 
second it is elliptical (EOA), both of them are not exact models, but EOA is 
very accurate for Earth effect prediction purposes. 

The systematic errors in the Sun position produced by the COA and the 
EOA affect (0) and (^ in two ways: i) they change the time in which the Sun 
cross the boundaries of a given Nadir angle interval, ii) they change h{t). The 
first way is the less relevant, since it mainly affects the residence time, which is 
the difference between two next crossing times, so that errors compensate. The 
second way is more relevant, because it affects the length of the neutrino path 
in a given geological structure. The analysis shows that EOA introduces errors 
in the Sun position smaller than 0.01°, while COA may introduce errors as 
large as few 0.01°. Most of these errors are periodic in nature, and they should 
average out. However, since Pg2 is a function of the solar position and the 
use of different sampling schemes (as Core, Mantle, Night, etc.) are equivalent 
to averaging over specific fractions of year only, it is not possible to test such 
approximations without numerical experiments, the most relevant case being 
the interplay between the change in flux due to the time dependence in the Sun 
- Earth distance R{t) and the time dependence in the Earth orbital velocity []T|. 
In conclusion, all the results presented here are obtained using the EOA, which 
assures an accuracy for < Pg2 > largely better than 1% for each kind of sample 



and in what regards the solar motion. 



3 Results 

At the basis of each time averaged probabihty < Pe2 > is the instantaneous 
one, Pe2- Figure 1 shows an example of the probability Pe2 for sin^ 2^^ = 0.01 
as a function of h and of the Resonance Density in the Earth pji (gr/cm^) []T|: 
pnYe = 6.57 X 10^{Am'^{ eV^)/E^{ MeV))cos2^y. There are two main peaks 
associated to resonance transitions in the core, C, and in the mantle M 
In addition, the spherical symmetry of the Earth implies that solar neutrinos 
will cross the resonance region twice. Since the neutrino oscillation length in 
matter, for solar neutrinos which undergo resonance transitions in the Earth, 
is of the order of the Earth radius, interference terms in Pe2 lead to oscillatory 
dependence of Pe2 on h for a fixed pr [Tg] . 

Increasing sm^26v, for sm^26v < 0.1, the qualitative properties of this de- 
pendences do not change significantly. The two peaks scales accordingly to 
sin^ 26v, with C increasing faster than M reaching a maximum value and be- 
gins to decrease while M is still increasing. As a consequence, at values of 
sin^2^y <0.13 most of the D-N asymmetry is generated by the MSW effect 
in the core. At larger sin^ 26v values the behaviour is more complicate, new 
peaks appears and (for sin^ 26v > 0.3) the bulk of the Earth transition is in the 
mantle. 

As Fig. 2a illustrates, the separation in Core and Mantle samples is very 
effective in enhancing < Pe2 > . The enhancement is particularly large for small 
mixing angles for which < Pe2 >^ can be up to six times bigger than < Pe2 >^ 
(for sin^ 2^y > 0.3 the core - enhancemente is only 30% -!- 40%). This suggests 
the possibility of a corresponding enhancement of the e~-spectrum distortions 
as well as of the energy integrated event rate. To have an enhancement in the 
averaged probability < Pe2 > is not enough to ensure an enhancement in the 
D-N asymmetry, especially at small mixing angles. The upper "window" of Fig. 
2b is a combined plot of Pq (dotted line), Pe2 (dashed line) and P© (solid line) 
for Core. In the lower "windows" there is a plot of the ratio: 

AAE.IAm ) = 2 p^(£^/^„2) ^ p^(£^/A„.) (4) 

which is the D-N asymmetry for the total survival probability. The large peak 
on the left side of the figure is an artifact of the presence of the adiabatic 
minimum in the probability P©. The bulk of the D-N asymmetry comes from 
the second peak in the figure. 



Figure 2 is an example of how each probabihty contributes in and in 
the related asymmetry. One of the important features that has to be taken 
into account in the discussion of is the position of the peaks in < Pp2 > 
with respect to the value of E^j l\vr? at which P© = 0.5 because it controls 
the sign and the magnitude of the D-N effect. It follows from eq. (|T]) that 
the asymmetry is zero each time P© = 0.5 irrespective of the value of < Pe2 >• 
If < Pe2 > > sin^^y, and P© < 0.5 (P© > 0.5), the asymmetry is positive 
(negative) which means the Earth effect increases [reduces) the night event rate 
The presence of a zero Earth effect can reduce the magnitude of the D-N 
asymmetry. The negative D-N effect is relevant for sin^ 2Qy < 0.004 0.006. 

The definitions for the recoil-e~ spectrum with (without) MSW effect S^{Te) 
{So{Te)), Tg being the recoil-e"kinetic energy and s = D, N, C, M; and 
for the energy integrated event rates with (without) MSW effect 7l^{Te^th) 
{lZQ{Tg^th)) , Tg^t/i being the recoil-e^threshold kinetic energy, in terms of the 
survival probability P©, are reported in [0| Here it is enough to recall the 
definition of the distortion for the recoil-e~spectrum due to the MSW effect: 
6S'{Te) = S'{Te)/So{Te), s = D, N, C, M, and the definitions for the D-N 
asymmetries for the recoil-e~spectrum ^|)_^(re) and the energy integrated 
event rate A^jj_^{Te^th)- 



with s = D, N, C, M. Both these observable quantities are solar model 
independent, and are related to A%{Ejy/ Am^), s = D, N, C, M, DC. For fixed 
values of sin^ 26v and Tg th'- 



N) 



a -[3 max(^^(P/Am^) + Pq^^^,^ 



(6) 



A^ (E-0 5MeV) AUE/Am')P^{E/Am^ 
Aj,_AE 0.5MeVj ~ a - (3 A^E / Am'') + Pr^E / Amy 



where s = D, N, C, M, DC, a = 0.190476, /3 = 0.0952381. In the upper for- 
mula ^|)_7v is regarded as a function of Am^, Pq is the value of Pq{E^/ Am'') 
computed for the E^/ Am' value which corresponds to the A^p{Ejj/ Am') max- 
imum max(^p). In the lower formula Am' is fixed, while it is assumed: 
E = Ey = Te + 0.5 MeV. Equations (0) overestimate the numerical results 
by a factor between 3% and 40%. The error decreases increasing sin^ 26v- 

One example of the predicted recoil - e^ spectrum distortion in the Sun, 
and one example of D-N asymmetries in the spectrum are shown in Figs. 3 



(see [0, |T3j for a larger set of examples) both the figures being computed for 
Ye = 0.467. The enhancement of the spectrum distortions and of the D-N 
asymmetry in the Core sample is clearly seen, especially for sin^ 2^y<0.013. 
The differences between the spectra computed for Y^, = 0.467 and for Y^ = 0.5, 
are largely negligible for the expected experimental accuracy. 

The differences between the Day and the Night {Mantle) sample - spectra 
are hardly observable, but for sin^ 29v > 0.008 and many values of Am^ from 
the "conservative" MSW solution region, ^^.^(Tg) may reach a value of 10% 
and even larger (20% — 30%). So, its measurement can further constrain Am^ 
and sin^ 29v values. 

Figures 4, 5 represent the iso-(D-N) asymmetry contour plots for the Night, 
Core and Mantle samples. Shown are also the "conservative" regions of the 
two MSW solutions of the solar neutrino problem (dashed lines) as well as the 
solution regions obtained in the reference solar model (inner thick dashed 
lines). Also the iso-(D-N) asymmetry contours for the Night and Core samples 
shown in Figs. 4 and 5 have been obtained both for Y^ = 0.467 (solid lines) and 
for Ye = 0.5 (dash-dotted hues). For given Am^ and sm'^29v the differences 
between A^_]si or Af^_j^ and A%_j^ are larger in the NA region. The magnitude 
of the D-N asymmetry in the AD region does not change significantly with the 
change of the sample. In contrast, in the NA region the D-N asymmetry in the 
Core sample is enhanced by a factor of up to six. The reduction in statistics 
for this sample reduces it to an effective 2.3 enhancement. From Fig. 5 it 
comes that the Core selection is a very effective method to enhance the D-N 
asymmetry despite the loss of statistics. 

The difference between the values of the Core sample asymmetry ^£_7v 
calculated for Y^ = 0.467 and for Ye = 0.5 depends strongly on the values of 
Am^ and sin^ 26v- However, Figs. 5 show that the main results illustrated here 
are stable against the Earth model uncertainty [0] . 

Because of the large differences between the core and the mantle structures, 
the Mantle and Core subsamples provide two independent measurements of the 
D-N effect. These can be combined to constrain better the neutrino parameters 
Am^ and sin^ 29v', utilizing the full statistics of the Night sample, especially in 
N A region . 

Figs. 5 suggests also the possibility to use the D-N effect to limit the 
neutrinos total flux []2j]. 



4 Conclusions 



In this talk the possibihty of a core-enhancement for the D-N effect in the Super 
- Kamiokande detector, proposed by example in JT^, was discussed in detail. 
It was remarked the need of a precise apparent solar motion reconstruction 
to produce the accurate D-N effect predictions required to study this subject. 
A computer code which fulfills this requirement was presented. Predictions for 
Super - Kamiokande detector observable quantities were obtained. It was shown 
that: I). The Core sample D-N asymmetry can be enhanced up to factor six 
compared to the Night and Mantle D-N asymmetries. The enhancement reduces 
to a factor 2.3 when the loss in statistics is properly accounted for. If Super - 
Kamiokande will be able to reach an accuracy of some percent in measuring the 
Core D-N asymmetry it is reasonable to expect it will be able to test this region 
through the D-N effect. II). The enhancement may be relevant for the recoil-e~ 
spectra distortion and the recoil-e^ spectra D-N asymmetry. III). The detection 
(or lack of detection) of a small D-N asymmetry may constraint the total 
neutrino flux from the Sun in a solar model independent way (presently this 
flux is uncertain for a factor two). IV). These results are qualitatively robust 
against the Earth model uncertainties. 



Acknowledgments 

The author is deeply indebted with Prof. W. Hampel and the organizing com- 
mittee of the Fourth International Solar Neutrino Conference for the kind con- 
cession of a financial support to attend to the conference. This research was 
partially supported by S.I.S.S.A. Trieste, Italy, research contracts: 4060 and 
5566. 



References 

[1] Q.Y. Liu, M. Maris and S.T. Petcov, Preprint SISSA 16/97/EP, January 
1997, E-archive report |hep-ph/9YU:^at)i| . 

[2] M. Maris and S.T. Petcov, SISSA Report ref. SISSA 17/97/EP, May 1997, 
E-archive report |hep-ph/'9V0b392| . 

[3] S.P. Mikheyev and A.Yu. Smirnov, Sov. J. Nucl. Phys. 6, 913 (1985); L. 
Wolfenstein, Phy. Rev. D 20, 2369 (1978). 

[4] J.N. BahcaU and M.H. Pinsonneault, Rev. Mod. Phys. 67, 781 (1995). 



[5] P.I. Krastev and A.Yu. Smirnov, Phys. Lett. B 338, 282 (1994); V. 
Berezinsky. G. Fiorentini and M. Lissia, ibid. B341, 38 (1994); N. Hata 
and P. Langacker, Phys. Rev. D52, 420 (1995). 

[6] P.I. Krastev and S.T. Petcov, reported by S.T. Petcov at the "Neu- 
trino'96" International Conference on Neutrino Physics and Astrophysics, 
June 13-19, 1996, Helsinki, Finland. 

[7] F.D. Stacey, Physics of the Earth, 2^'^ edition, John Wiley and Sons, 
London, New York, 1977. 

[8] J.N. Bahcall et al., Phys. Rev. D 54, 411 (1996). J.N. BahcaU, M. 
Kamionkowski and A. Sirlin, Phys. Rev. D 51, 6146 (1995). 

[9] S.P. Mikheyev and A.Yu. Smirnov, in Massive Neutrinos in Astrophysics 
and in Particle Physics, Proceedings of the Moriond Workshop, Tignes, 
Savoie, France, 1986, edited by O. Fackler and J. Tran Thanh Van (Edi- 
tions Frontieres, Gif-sur-Yvette, France, 1986), p. 355; in New and Exotic 
Phenomena, Les Arc, Savoie, France, 1987, edited by O. Fackler and J. 
Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, France, 1987), p. 
405. 

[10] A.J. Baltz and J. Weneser, Phys. Rev. D 50, 5971 (1994). 

[11] S.T. Petcov, Phys. Lett. B 200, 373 (1988). P.I. Krastev and S.T. Petcov, 
Phys. Lett. B 205, 84 (1988). 

[12] B. Bertotti and M. Maris, Variabilitd Temporale del Flusso di Neutrini 

Solari in Esperimenti Sotterranei (Neutrino Flux Time Variability in Underground 

Experiments), reported in the Ph. D. thesis of M. Maris, Department of The- 
oretical and Nuclear Physics (D.F.N.T.), Pavia University, 1996, Italy. 

[13] M. Maris and S.T. Petcov, E-archive report [hep-ph/9YU51^UY| , February 
1997. 



Pe2 in the Earth 




D-N Asymmetry: Night Sample 




Figure 4 



D-N Asymmetry 




sin 26 V 



Figure 5 



sin'(20j = 0.0100 



0.4 



0.3 



0.2 



<Pe2> 



Night, Core, Mantle 
Deep Core 




p, (gr/cm ) 




0.6 
0.4 
0.2 



si: 



1 O^To^To^ryTo^K^K) ^ ° 
E./Am^ (MeV/eV^) 



Figure 2. a) < Pe2 > for Night (short dashed hne), Mantle (dotted hne), Core (black hne) 
and Deep - Core (long dashed line) samples, b): up frame (black line), < Pe2 > (dotted 
line) and (dashed line) for the Core sample; down frame Ap{E^ / dms) for the Core 
sample. 
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Figure 3. Recoil - e~ spectrum distortion (upper) and D-N asymmetry (lower), for Day (long 
dashed line). Night (short dashed line). Mantle (dotted line). Core (black line) samples. 



